Controlling soliton interactions in Bose-Einstein condensates by synchronizing the 
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We present how to control interactions between solitons, either bright or dark, in Bose-Einstein 
condensates by synchronizing Feshbach resonance and harmonic trap. Our results show that as long 
as the scattering length is to be modulated in time via a changing magnetic field near the Feshbach 
resonance, and the harmonic trapping frequencies are also modulated in time, exact solutions of 
the one- dimensional nonlinear Schrodinger equation can be found in a general closed form, and 
interactions between two solitons are modulated in detail in currently experimental conditions. 
We also propose experimental protocols to observe the phenomena such as fusion, fission, warp, 
oscillation, elastic collision in future experiments. 
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I. INTRODUCTION 

Observation of Bose-Einstein condensates (BECs) in 
gases of weakly interacting alkali-metal atoms has stim- 
ulated intensive studies of the nonlinear matter waves. 
One of the central questions in this field is how to explore 
properties of BECs. It is known that interatomic inter- 
actions greatly affect a number of properties of BECs, in- 
cluding both static (such as the size, shape and stability) 
and dynamic ones (like the collective excitation, soliton, 
and vortex behavior, etc.). A common practice to change 
the interaction strength, even its sign, is to modulate the 
s-wave scattering length, a s , by using the Feshbach reso- 
nance with a tunable time-dependent magnetic field B(t) 
[E 1 1 11 0, 1 S |: a s (t) = a 00 [l- A/(B(t) - B )}, 
where is the off-resonance scattering length, t is the 
time, Bq and A are the resonance position and width, 
respectively. This offers a good opportunity for manipu- 
lation of atomic matter waves and nonlinear excitations 
in BECs. In real experiments, various forms of the time 
dependence of B(t) have been explored [1, lift [IE ' obser- 
vation of dark and bright solitons have been reported. In 
theoretical studies, several forms of time- varying scatter- 
ing lengths have been proposed and treated separately, 



such as the exponential function e A * Tp, 13 L 14 1, or the 
periodic function g + gi sin(fit) [H, lH, and 



bitrary time- varying scattering length a s (t) , then further 
discuss how to control the interaction of solitons includ- 
ing the bright and dark solitons. We observe several in- 
teresting phenomena such as fusion, fission, warp, oscil- 
lation, elastic collision in BECs with different kinds of 
scattering length correspond to different real experimen- 
tal cases. 



II. THE MODEL AND SOLITON SOLUTIONS 



Consider condensates in a harmonic trap V(r) = 
mut^y 2 + z 2 )/2 + muj 2 x 2 /2, where m is atomic mass, oj± 
and iO\ the transversal and axial frequency, respectively. 
Such a trap can be realized, for instance, as a dipole trap 
formed by a strong off-resonant laser field. In the mean- 
field theory, the dynamics of BEC at low temperature 
is governed by the so-called Gross-Pitaevskii (GP) equa- 
tion in three-dimensions. If u>± !3> |^i|, it is reasonable 
to reduce the GP equation for the condensate wave func- 
tion to the quasi-one-dimensional (qua si-lD) nonlinear 
Schrodinger equation [U HE H US Q , 



.dih 



ld 2 ip 
Ydx 2 



a B 



x ip, 



(1) 



In the present paper, we will consider the general case 
with arbitrary, time- varying scattering length a s (t), and 
discuss how to control dynamics of solitons in BECs 
by synchronizing the Feshbach resonance and harmonic 
trap in current experimental conditions. We first ob- 
tain a family of exact solutions to the general nonlinear 
Schrodinger equation with an external potential and ar- 



where the time t and coordinate x are measured, respec- 
tively, in units of ujJ 1 and aj_, with a± = ^h/muj±; iff 
is measured in units of l/{y/2 / Ka 2 1 a,B), with as as the 
Bohr radius. The key observation of the present paper is 
that if we allow the axial frequency of the harmonic trap 
becomes also time dependent u>\ — u>i(t), and require 
it to satisfy the following integrability relation with the 
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scattering length a s (t), we have 

1 d 2 a s {t) 2 da s (t) u 2 (t) _ 
a s (t) dt 2 a 2 (t) [ dt ' u>l ' { ' 

Then, the nonlinear GP equation {T]) with time- varying 
coefficients can be reduced to the standard nonlinear 
Schrodinger equation and exactly solved, with the fol- 
lowing general solution: 

t 

ip(x, t) = exp[- J r(t)dt]<f>(X, T) exp\iT(t)x 2 }, (3) 

Co 

where <f> is an arbitrary function of X and T, 

with new spatial and temporal variables X — 

t 2 t t 

Aiexp[-2 fT(t)dt]x, T = 4f- / exp[~4 J T(t')dt']dt. A 1 

to to to 

is a real constant, which together with T(t) are deter- 
mined by 

t 

a s (t) = ±a B A 2 1 exp[-2 J T{t)dt}. (4) 

to 

Meanwhile, the trapping frequency can also be ex- 
pressed in terms of T(t) given by 

^M = -2iwr(t) 2 . (5) 

From Eqs. (0| and ©, we see that both the scattering 
length and the trapping frequency can be expressed in 
terms of the T function. It shows that the trapping po- 
tential can become repulsive, we will give a detailed expo- 
sition of this situation in the following section. Since the 
trapping potential we consider here is the cigar-shaped 
harmonic potential (here-after, the frequency co± is not 
varying), once the function of scattering length is de- 
termined, the r function and the function of the trap 
potential can also be determined. Note that the exact 
solution tj}(x, t) can be obtained for arbitrary time de- 
pendence of a s (t) , since we can always choose an appro- 
priate time-dependent axial frequency, wi(t), to satisfy 
the integrability relation. 

When the interatomic interaction is attractive, i.e., 
a s (t) < 0, Eq. (P) has bright iV-soliton solutions. The 
simplest case for studying soliton interactions is the two 
bright solitons solution, for which ip(x,t) is expressed by 
Eq. © with <j)(X, T) given by 

<j>{X,T) = 2bcxp{i[cX - (c 2 - b 2 )T- ip]}x 

(2b 2 T-i) cosh[b(X-2cT-X n )]+ib(X-2cT) sinh[b(X -2cT-X )] 
costi A [b(X-2cT-X n )] + b 2 [(X-2cT) 2 +4b 2 T] ' 

; (6) 

where 6, c, tp, Xq are arbitrary constants. 

When the interatomic interaction is repulsive, i.e., 
a s (t) > 0, there are dark JV-soliton solutions to Eq. fTJ). 



For the two dark solitons solution, ip(x,t) is again ex- 
pressed by Eq. ([3]), but with <j>(X,T) given by 

<t>(X, T) = ri exp{i[hX - (I 2 - 2\ Tl \ 2 )T + 5t}} x 

l+e[Zi cx.p(£ 1 )+Z 2 cxp(g 2 )]+e 2 Ai 2 ZiZ 2 cxp(£i+g 2 ) ( 7 \ 
l+£[c X p(£ 1 )+cxp(£ 2 )]+^A 12 c X p(C 1 +? 2 ) ' \') 

where Zt = (-fii+V 4 l r i| 2 - ^ 2 )/(^i+ V 4 l r i I 2 = 

z 2 = (-o 2 + V 4 lnl 2 - n 3)/( fi 2 + y4|riP- n|), 
A ig [(sit - si 2 ) 2 + {^ A\ Tl \ 2 ~Sl 2 - 

y g^F^of ) 2 ]/^ + Sl 2 ) 2 + {^ A\ Tl \ 2 -sij - 
y/^-SXif], jt = WM t^-SI 2 - 2l 1 )Q 1 T + 
SliX + t[°\ & = (y/4\ Tl \ 2 -SI 2 - 2h)Sl 2 T + Sl 2 X + £_ { 2 °\ 

It, St, Sit, Sl 2 , \ S} 2 \ e are real constants, and t\ is 
complex constants. 



III. EFFECTS OF THE TIME-DEPENDENT 
MAGNETIC FIELDS ON THE SOLITONS 

Now we consider the elementary applications of solu- 
tions ([3]) with ((6|) and ([7]) respectively, with linear, expo- 
nential and sinusoidal time dependence of the magnetic 
field via Feshbach resonance, and propose how to control 
dynamics of solitons in BECs by synchronizing Feshbach 
resonance and harmonic trap in future experiments. 

A. MAGNETIC FIELD RAMPED LINEARLY 
WITH TIME 

In real experiments [H, [Io| . the magnetic field is lin- 
early ramped down with time t. We can design an ex- 
perimental protocol to control the soliton interaction in 
BECs near Feshbach resonance with the following steps: 
(i) In real experiment of 87 Rb atoms, the scattering 
length can be chosen as a function of magnetic field, i.e., 
a s (t) = aoo[l — A/(B(t) — Bo)] Q, where the off-resonant 
scattering length Ooo = 108as, as is the Bohr radius, Bq 
is the Feshbach resonance position, and A is the res- 
onance width, respectively. The best-fit value for the 
width is A = 0.20 G, resulting in B = 1007.40 G. Near 
the Feshbach resonance, the field B(t) varies linearly with 
the rate 0.02 G/ms. For a better understanding, we plot 
Fig. 1, which shows the scattering length and the trap- 
ping frequency vary with time when the field B(t) ap- 
proaches the Feshbach resonance position Bo. (ii) The 
realistic experimental parameters for a quasi-lD repul- 
sive condensate can be chosen N ~ 10 3 atoms and with 
peak atomic density no = 10 5 c7to _1 . Then the scattering 
length a s is of order of nanometer, e.g., \a s \ = 5.8 nm 
for a 87 Rb condensate, and uj± = 2tt x 400 Hz, with the 
ratio u>t/u>x being very close to zero. 

The validity of the GP equation relies on the condi- 
tion that the system be dilute and weakly interacting: 
n|a s | 3 -C 1, where n is the average density of the conden- 
sate. Applying the above conclusions to real experiments, 
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FIG. 1: (color online) The scattering length (red dotted line) 
and the trapping potential (blue solid line) vary with time 
when the magnetic field ramps linearly with time t. (a) The 
selected field B(t) varies from 1007.54 G to 1007.42 G, corre- 
sponding to the attractive interaction between atoms, (b) The 
selected field B(t) varies from 1007.70 G to 1008.50 G, cor- 
responding to the repulsive interaction between atoms. The 
scattering length is measured in units of Ooo, the axial fre- 
quency is measured in units of w\, and the field B(t) varies 
linearly with the rate 0.02 G/ms. 



we need to examine whether the validity condition for the 
GP equation can be satisfied or not. In the ground state 
for 87 Rb condensate, the scattering length is known to be 
\a s \ = 5.8 nm [§]; the typical value of the density ranges 
from 10 13 -10 15 cm" 3 . So n\a s \ 3 < 1CT 3 < 1 is satisfied. 
Moreover, the experimental data agree reasonably well 
with the mean-field results 25[ , which further proves the 
validity of the GP equation with \a s \ = 5.8nm. Another 
important issue is quantum depletion of the condensate, 
which is ignored in the derivation of the GP equation. 
The physics beyond the GP equation should also be very 
rich, and we will work on more rigorous solutions beyond 
the GP equation in the future. 

Figure 2(a) shows the bright solitons interaction in 
BECs near the Feshbach resonance. In this case, the 
s-wave scattering length a s < 0. According to the in- 
tegrability relation, the time-dependent axial frequency 
is imaginary which indicates a repulsive potential. As 
shown in Fig. 1(a), when the field B(t) approaches the 
Feshbach resonance position B , the absolute value of 
the scattering length increases, but the time-dependent 
axial frequency decrease linearly. With the increasing 
of the absolute of the scattering length, the interactions 
between atoms become stronger, the peak of each soli- 
ton increases and its width decreases. Meanwhile, under 
the expulsive potential, the two bright solitons will be 
set into motion. As a result, the left bright soliton feels 
two forces which come from the right bright soliton and 
the trapping potential, it drives the left bright soliton to 
the right-hand side. The right-hand one moves slower 
than the left-hand one. Finally, the distance between 
the solitons becomes smaller. When the field B(t) in- 
finitely approaches Bo, the absolute value of the scatter- 
ing length and axial frequency become infinite. The two 
solitons interact very strongly and almost merge into one 
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FIG. 2: (color online) Controlling matter wave bright and 
dark soliton interaction near the Feshbach resonance (Bo = 
1007.40 G) when the magnetic field ramps linearly with time 
t. (a) The selected field B(t) varies from 1007.54 G to 1007.42 
G. The parameters are given as follows: b = 0.01, Xo = 3000, 
c = 0.002 in Eq. © with Eq. ©. (b) The selected field B(t) 
varies from 1007.70 G to 1008.50 G. The parameters are given 
as follows: |n| 2 = 0.0001, fii = 0.01, fi 2 = -0.02, £< 0) = 112, 
0, e = 1 in Eq. © with Eq. ©. 
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with a very high peak and the narrowest width. After at 
least close to 1007.50 G, the absolute value of the atomic 
scattering length becomes | a s (t) | = 5.7 nm < 5.8 nm for 
quasi-lD 87 Rb gas mentioned above. This means that the 
stability of soliton and the validity of ID approximation 
is maintained from 1007.54 G to 1007.50 G. With further 
increasing of \a s \ [for example, while B(t) — 1007. 46G, 
\a s \ should be 13.3nm], the system may be beyond the va- 
lidity of the GP equation. Therefore, the phenomena dis- 
cussed in Fig. 2(a) should be observable within the cur- 
rent experimental condition from 1007.54 G to 1007.50 
G. With synchronized Feshbach resonance and harmonic 
trap to change the scattering length and axial frequency, 
we can easily control matter wave soliton interactions 
and obtain a new type of atom laser with manipulatable 
intensity. 

Figure 2(b) shows the interaction of dark solitons with 
magnetic field being selected in the range from 1007.70 G 
to 1008.50 G. In this shown in Fig. 1(b), the s- 

wave scattering length a s > 0, and is in proportion to the 
magnetic field, the axial time-dependent decrease linearly 
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and asymptotic approaches to zero, it means that the sys- 
tem becomes a self-confined condensate. Initially, there 
is only one dark soliton in BECs. Repulsive interaction 
between atoms becomes stronger when the absolute value 
of a s increases. This causes the dark soliton to split and 
finally become two dark solitons, meanwhile, each soliton 
increases its peak and compresses its width. However, 
contrary to the bright soliton shown in Fig. 2(a), one 
dark soliton moves away from the other. After at least 
close to 1008.50 G, the absolute value of the atomic scat- 
tering length becomes |a s (i)| = 4.67 nm, which is safely 
smaller than 5.8 nm for quasi-lD 87 Rb condensate, thus 
the stability of the soliton and validity of the ID ap- 
proximation are maintained. We conclude that the dark 
soliton fission phenomenon revealed here can be realized 
under the current experimental condition. 



B. MAGNETIC FIELD VARYING 
EXPONENTIALLY WITH TIME 

When the field is exponentially ramped down as 
exp(— t/r) to a selected field between 545 G and 630 G 
Q, where r = 40 ms, a s small and negative or small 
and positive, the interaction parameter g(t) near the res- 
onance varies exponentially with time: g(t) = a s /as = 
±0.01 exp(At), where A = |k>i|/wj_ <C 1. The integrabil- 

2 

ity relation reads —A 2 = — 1 », and can be satisfied auto- 
matically as long as the time-dependent axial frequency 
is imaginary which indicates a repulsive trapping poten- 
tial. With the same parameters as in the experiment [26j ] . 
i.e., N ps 10 3 , u± = 2?r x 700 Hz, wi = 2iiT x 21 Hz, and 
A = 0.03, the interactions between the two bright soli- 
tons is shown in Fig. 3. It is interesting to observe that 
in the expulsive parabolic potential, the bright solitons 
are set into motion and propagate in the axial direction. 
With time going on, |a«(i)| increases. We could observe 
an increase in their peaking values and a compression in 
widths, besides that the spacing between them decreases 
after they are initially generated on different positions in 
the trap, which is evidence for a short-range attractive 
interaction between solitons. Finally, they almost merge 
and fusion. This phenomena is different from the case 
of [^, [2f|. There, the bright solitons are set in motion 
by off setting the optical potential and propagate in the 
potential for many oscillatory cycles with the period 310 
ms, the spacing between the solitons increase near the 
center of oscillation and bunches at the end points. The 
difference is mainly caused by two factors: one is the 
time-varying scattering length strongly affect the inter- 
action between the solitons, and the other reason is the 
repulsive force provided by the potential. Meanwhile, 
with the increasing of the scattering length, the attrac- 
tive interaction between atoms become stronger, this will 
leads to an attractive interaction between solitons. Even- 
tually, the outcome of these two forces will determine the 
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FIG. 3: (color online) Controlling matter wave bright soliton 
interaction when field varying exponentially with time. The 
parameters are given as follows: A = 0.03, b = 0.5, Xo — 20, 
c = 7 in Eq. (|3} with Eq. (|6]). The color corresponds to the 
BEC density, with yellow (gray) being the smallest and blue 
(dark) being the largest. 



motion of the two bright solitons. 

Fusion is very interesting phenomenon and it comes 
from the interatomic attractive interaction. In other 
words, with time going on, both bright solitons change 
their positions, warp in a certain radian, and almost 
merge into one single soliton. Such morphology has been 
observed in coronal plasma 27j . We hope that such mor- 
phology would be detected in BEC experiments too in the 
near future. 

In real experiments , the length of the background of 
BECs can reach at least 2L = 370 fim. At the same time, 
in Fig. 3, solitons travel from x — 50 to 200, i.e., 150 x 1.4 
fim= 210 fjjxi. [The dimensionless unit of the coordinate, 
Ax — 1, corresponds to a±_ — (h/muj^) 1 ^ 2 =1.4 /zm]. 
We indeed have 210 /im < 370 /im, a necessary condi- 
tion for observing the morphology in BEC experiments 
[H |. Additionally, after at least up to 100 dimensionless 
units of time, |a s (t)| reaches the value 0.2as, which is 
less than |a/i no i| = 4as. This means that during the 
time evolution, the stability of solitons and the validity 
of ID approximation can be maintained as displayed in 
Fig. 3. Therefore, the phenomena discussed in this case 
are also expected to be observable within the current ex- 
perimental capability. 

The interactions between two dark solitons are also in- 
triguing. The first experimental evidence of attraction 
between dark solitons in nonlocal nonlinear media has 
been presented [28]] . Our results (3), (5), and (7) also 
indicated that attraction between dark solitons should 
be observable in BECs with repulsive long-range inter- 
atomic interaction. In a previous case, the field is lin- 
early ramped down as time, leading to the trap axial fre- 
quency u>i = and external potential vanishing. It means 
that the system becomes a self-confined condensate. So 
the repulsion between the dark solitons is observed in 
BEC with repulsive interatomic interaction [seen in Fig. 
2(b)]. However, in the present case, the field varies ex- 



5 



1.5i . . . . 

VWWWWWVW\A 

0.5 ■ 

a (t) m=0.1 

s o 

■ ■ ■ m=0.01 

-0.5 m=0.3 
-1 

-1.5I 1 1 1 1 1 

10 20 t 30 40 

FIG. 4: (color online) The scattering length varies with time 
when field varying periodically with time t. The yellow curve 
(gray) denotes the attractive interaction case, at t — 30 
(marked by an arrow), the amplitude of the scattering length 
is larger than before by suddenly jumping the amplitude of 
the ac drive mtom = 0.3, w — 30. The red curve denotes 
the repulsive interaction for m = 0.1 (dark solid line) and 
m = 0.01 (dark dotted line), w = 2. All of the lengths are 
measured in units of as. 

ponentially with time, the trap axial potential is not 
vanishing and time independent, the system is in an ex- 
pulsive parabolic potential. Following the experimental 
setup in [20|, we can first create a BEC in the quasi- 
one-dimensional potential. Second, the trap potential is 
tuned to the value in our paper, meanwhile, the scat- 
tering length varies exponentially with time to a small 
and positive value. The potential provides an attractive 
force which counters the natural repulsion of the solitons. 
Finally, the competition between these two forces will de- 
termine the outcome of the interaction between solitons. 
If the attractive force which is caused by the potential is 
stronger than the natural repulsion of the solitons, the 
two dark solitons will move towards the other, which is 
the evidence for a short-range attractive interaction be- 
tween dark solitons. However, the lifetime of a BEC in 
current experiments is of the order of 1 s and the region 
of a BEC is small, the solitons will be dissipation in their 
motion before reaching the edge of the condensate. 

C. MAGNETIC FIELD VARYING 
PERIODICALLY WITH TIME 

It was observed that a small sinusoidal modulation 
of the magnetic field close to the Feshbach resonance 
gave rise to a modulation of the interaction strength, 
g(t) = a s (t)/ciB = ±[l+msin(orf)] fill ], where the ampli- 
tude m of the ac drive was small and satisfied < m < 1. 
According to the integrability condition, the axial fre- 
quency of the harmonic potential should be 

2 2 

w iW = muJ U f [s\n(ut)+m+mcos 2 {ut)]. (8) 

[1 + m sm(ujt)\ z 

Now, we investigate how the amplitude of the ac drive 
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FIG. 5: (color online) The trapping frequency varies with 
time when field varying periodically with time t. (a) m — 0.01 
(0 < t < 30), m = 0.3 (30 < t < 50), and w = 30 for the case 
with atomic attractive interaction, (b) m = 0.01 (red solid 
line), m = 0.1 (blue dashed line), and w = 2 for the case with 
atomic repulsive interaction. The frequency is measured in 
units of Wj_. 

can be used to control the bright soliton interactions. 
As shown in Fig. 6, for the case with atomic attrac- 
tive interaction, when the amplitude is small, m = 0.01 
(0 < t < 30). The periodically varying of the scattering 
length and the trapping frequency are small as shown in 
Figs. 4 and 5, the repulsive and attractive force between 
solitons can become balanced. Two bright solitons will 
move in parallel with their separation keeping constant. 
This property would be interesting for optical communi- 
cation with low bit-error rates. In BEC, this effect may 
play an important role in potential application of matter 
wave communication with atom lasers. If the amplitude 
is increased to the value m = 0.3 (30 < t < 50), the 
solitons begin to oscillate due to the temporal periodic 
modulation of the s-wave scattering and trapping poten- 
tial are both stronger than before. This phenomena is 
very similar to the evolution of optical solitons. 

We also studied the effect of ac drive on interactions 
between the dark solitons. First, the amplitude m of the 
ac drive is chosen small, m = 0.01, this leads to a small 
periodical modulation of the scattering length and the 
trapping frequency. The repulsive interaction between 
atoms mainly leads to the formation of the dark solitons, 
this force cannot lead to a oscillation of the dark soli- 
tons due to the small change in the scattering length. As 
shown in Fig. 7(a), a faster soliton is generated behind 
a slower one in the mutual moving direction, after an in- 
terval of time, the faster one pulls up to the slower one 
and their elastic collision happens. After the collision, 
the faster one passes through the slower one and their 
parameters did not change, which remarkably indicates 
no energy exchange between the two dark solitons. For 
larger ac drive, for example, m = 0.1, which is 10 times of 
the value of the previous m, with the same initial condi- 
tion, because of the stronger modulation of the scattering 
length and the trapping frequency, both solitons move 
forward, meanwhile, they will oscillate back and forth, 



FIG. 6: (color online) Controlling matter wave bright soliton 
interaction when field varying periodically with time, where 
(a) m = 0.01 (0 < t < 30), (b) m = 0.3 (30 < t < 50). The 
other parameters are as follows: u) = 30, 6 = 1, Xo = —9, 
c = 0.4 in Eq. © with Eq. ©. 

as shown in Fig. 7(b). In all of the above cases, the trap 
potential can become repulsive during the entire process, 
but the attractive potential is stronger than the repul- 
sive one during a period, the outcome of the potential is 
an attractive one which can be seen in Fig. 5. It shows 
that a BEC with repulsive interaction between atoms is 
confined in the trap, this is different from the case B. 

Another important problem is the atom loss. As the 
external magnetic field is driven close to the resonant 
value, the rate of loss of atoms is increasing rapidly in 
the vicinity of the Feshbach resonance, while only a small 
fraction of atoms remain as soliton. In all of the above 
cases, the scattering length is small, the validity of the 
GP equation is satisfied. Meanwhile, when the trapping 
potential is modulated according to the integrability re- 
lation, both the rate of the untrapped atom and the col- 
lective excitations will be further discussed. 
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FIG. 7: (color online) Controlling matter wave dark soliton 
interaction when field varying periodically with time, where 
(a) m = 0.01, (b) m = 0.1. The other parameters are as 
follows: u> = 2, |nj 2 = 1, fii = 1.2, Q 2 = -1.8, h = 0.02, 
= -0.01, Ca 0) = -0.2, e = 1 in Eq. © with Eq. ©.The 
red (with smaller slope) and blue correspond to the slower 
one and the faster one respectively. 
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IV. CONCLUSION 

In summary, we express how to control soliton inter- 
action in BECs with arbitrary time-varying scattering 
length in a synchronized time-dependent harmonic trap. 
When the integrability condition is satisfied, we obtained 
the exact solutions analytically, and explored the inter- 
action of the bright and dark solitons in BECs with Fesh- 
bach resonance magnetic field linearly, exponentially, and 
periodically dependent on time. In these typical exam- 
ples, we find several interesting phenomena involving soli- 
ton interactions, such as fusion, fission, warp, oscillation, 
elastic collision, etc. 2^, 30j]. We further discussed how 
to control interactions between bright or dark solitons, in 
BECs in realistic situations, which allows for experimen- 
tal test of our predictions in the future. These phenom- 
ena open possibilities for future applications in coherent 
atom optics, atom interferometry, and atom transport. 
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